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Abstract. We consider an implicit finite difference sclieme on uniform grids 
in time and space for tfie Caucliy problem for a second order parabolic stoclias- 
tic partial differential equation where the parabolicity condition is allowed to 
degenerate. Such equations arise in the nonlinear filtering theory of partially 
observable diffusion processes. We show that the convergence of the spatial 
approximation can be accelerated to an arbitrarily high order, under suitable 
regularity assumptions, by applying an extrapolation technique. 



1. Introduction 

Motivated by the nonhnear filtering theory, we consider the Cauchy problem for 
the stochastic partial differential equation (SPDE) 

di 

(1.1) du = {a^^D^Dpu + /) dt + J2 (b'^^D^u + gP) dwP 

p=i 

with initial condition u(0) — uq where (u''')^]^j^ is a di-dimensional Wiener mar- 
tingale for integer di > 1 and the summation convention is used with respect to 
the repeated indices a, 13 G {0, for integer d > 1. Here Da ■— d/dxa, for 

a G {1, . . . ,d}, denotes the spatial derivative with respect to the direction a and we 
let Do be the identity. A special case of this equation, when the free terms / and g 
vanish, arises as the equation for the unnormalized conditional density of a signal 
process with respect to an observation process in the nonlinear filtering theory and 
is often referred to as the Zakai equation (see [131 IIZl HI])- The behavior of this 
equation is governed by the quadratic form 

d 

a, 13 = 1 

for A"/^ := 2a"P - b°'Pb'^P and z eR''. In [H], it is emphasized that in the setting 
of the nonlinear filtering theory one is only guaranteed the nonnegative definiteness 
of the matrices A, that is, when (11.11) satisfies a degenerate stochastic parabolicity 
condition (c/. [TTJ [TB] where the solvability of this equation is studied under the 
uniform nondegeneracy of the matrices A) . In application, these problems are high 
dimensional in nature and the solutions are required on-line. Therefore accurate 
and efficient numerical methods are desired for solving the Cauchy problem for p.ip 
under a degenerate parabolicity assumption. 

The present manuscript concerns the accuracy of a space-time, that is, a fully 
discretized, finite difference scheme on uniform grids in time and space for ap- 
proximating the solution to the Cauchy problem for (|l.ip under the degenerate 
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parabolicity assumption. In general the rate of convergence of finite difference 
schemes is known to be unsatisfactory in high dimensional settings. We prove that 
the rate of convergence of the spatial approximation for our space-time scheme can 
be accelerated to an arbitrarily high order with respect to the computational effort 
by applying an extrapolation technique. That is, we show that the rate of the 
strong convergence of the spatial approximation to the temporal discretization can 
be accelerated to any order of accuracy if the initial conditions, coefficients, and 
free terms are sufficiently smooth in space and the matrices A can be decomposed 
as 

for matrices a sufficiently smooth in space. While the requirement that the A admit 
such a decomposition is quite restrictive, this condition is satisfied in the nonlinear 
filtering problem even in the general case of correlated signal and observation noises 
when the diffusion coefficients of the signal noise are sufficiently smooth. 

The extrapolation technique that we employ to obtain higher order convergence 
is often referred to as Richardson's method, after L.F. Richardson who suggested 
the method for accelerating the convergence of finite difference approximations for 
certain partial differential equations (PDEs) (see [IHIIIS]). The method relies on 
the existence of an asymptotic expansion for the error between the approximate 
and true solutions to a continuous problem in powers of the discretization param- 
eter. Richardson observed that by taking appropriate weighted averages of the 
approximation at different mesh sizes certain lower order terms in the expansion 
vanish yielding a higher order rate of convergence. Therefore, it is important to 
give sufficient conditions under which such expansions exists. We emphasize that 
not only does the existence of the asymptotic expansion allow us to apply Richard- 
son's method to an arbitrarily high order, but also it allows us to measure the rate 
of convergence in the supremum norm. Richardson's method has been thoroughly 
studied in the literature, see for example the book [14] which provides a study of 
Richardson's method for finite difference schemes for deterministic PDEs and the 
survey articles [H [TU] on convergence acceleration methods. Part I of the book [52] 
concerns Richardson's method and algorithms for its implementation for PDEs; our 
results are of a more theoretical nature and much work still needs to be done to 
implement them. 

While finite difference schemes for PDEs and, to a lesser extent, for SPDEs are 
well studied in the literature (for example, see [221 IS] and the references therein) 
there are only a few results for degenerate parabolic equations and even fewer 
results concerning convergence acceleration for degenerate equations. Sharp rates 
of convergence are given in [3] for monotone finite difference schemes for possibly 
degenerate parabolic and elliptic deterministic PDEs. In [6] Richardson's method 
is applied to monotone finite difference schemes for possibly degenerate parabolic 
deterministic PDEs to accelerate the rate of convergence. Recently, in [1] a rate of 
convergence is given for a class of finite difference methods, that approximate in 
space via finite differences while allowing time to vary continuously, for degenerate 
parabolic SPDEs and sufficient conditions are given for accelerating the rate of 
convergence for the approximation in space. 

The current manuscript extends the results of [1] to a fully discretized scheme. 
We also mention [5], where results similar to those of the present manuscript are 
given under the strong parabolicity condition. A principle contribution of the 
present work is to provide estimates in the supremum norm in appropriate spaces 
for the solutions to the space-time scheme and the discretization in time under the 
degenerate parabolicity condition. The methods used to provide the requisite esti- 
mates in [3] are not tenable in the discrete time case. Further, we mention that we 
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have chosen to consider here imphcit schemes as we beheve that these are favored 
from a practical standpoint and because such schemes are unconditionally stable. 
We note that, to the author's knowledge, there are no results that give the rate of 
convergence of the implicit time scheme for SPDEs under a degenerate parabolic- 
ity assumption and it will be the subject of a future work to give such a rate of 
convergence for the implicit time scheme as well as more general methods. 

The paper is outlined as follows. In the next section, we begin by presenting 
our time scheme and our space-time scheme for approximating the solution to the 
Cauchy problem for as well as some preliminaries and assumptions. We then 
state our main results. Theorem 12.131 gives sufficient conditions for the existence of 
an asymptotic expansion for the error between the space-time approximation and 
the temporal discretization in powers of the spatial mesh size. Theorem 12.141 gives 
sufficient conditions for a generalization of Theorem 12.131 namely, the existence of 
such an expansion for differences of the solution. Then Theorems 12.151 and 12.171 
using the aforementioned expansions, give an accelerated rate of convergence for the 
spatial approximation and for derivatives of the spatial approximation, respectively. 
The proof of Theorem 12.141 and hence Theorem I2.13[ is given in Section 0] after 
some preliminary estimates are proven in Section [31 

We end this section by introducing some notation that will be used through- 
out this work. For integer d > 1, let R'' be the space of Euclidean points x = 
(xi, . . . ^Xd)- We denote the cr-algebra of Boreal subsets of R*^ by B^R!^). Recall 
that we denote by Da '■— d/dxa for a € {1, . . . ,d} the spatial derivative with re- 
spect to the direction a and let Dq be the identity. For an integer to > 0, we denote 
by :— Vl^™(R'') the usual Hilbcrt-Sobolev spaces of function on R"*, defined as 
the closure of C^(R'') functions cj) : R'^ K'^ in the norm 




where D'' — D^^ . . . D^'' for a multiindex p = (pi, . . . , pd) of length \p\ — pi+- ■ ■+pd- 
For an integer s > 0, we will use the notation D'^cj) to denote the collection of all 
sth order spatial derivatives of (j), that is, D'^cf) := {DP<j) : \p\ < s} for functions 
(j) = (t){x) for X e R''. We note that for := L'^{W^) = W§ we will denote 
the norm by || • ||o and we will use (•, •) to denote the usual inner product in that 
space. Let {^l,J^,P) be a complete probability space and let J-{t), t > 0, be an 
increasing family of sub-cr-algebras of J- such that J-{0) is complete with respect 
to (J", P). For a fixed integer di > 1 and a constant T G (0, oo) let {'w'')^^^ be a 
given sequence of independent Wiener processes carried by the complete stochastic 
basis {^l, T, {T{t))t>o, P). For the fundamentals of the nonlinear filtering theory, 
we refer the reader to the book yy and for basic notions and results from the theory 
of SPDEs we refer the reader to [H] . 

We collect the following notation for our discretizations and differences. For 
fixed T £ (0,1), we partition [0,r] into a fixed integer n > with mesh size r 
obtaining the time grid 

{it; i G {0,1, . . . , n}, rn — T} . 

We define (jji :— (j){iT) for functions <f> depending on i e [0,T]. In particular for 
i € {I, . . . ,n}, we will use the notation 

to denote the increments of the Weiner process for each p E {1, . . . ,di} and J^i 
J^(ir) to denote the filtration. For h E R\{0} and a finite subset A C R'' containing 
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the origin we define the space grids 

Gh := {Xih + • • • + \ph-p e {1, 2, . . . }, A, e A U (-A)} 

with mesh size \h\. We denote Aq := A \ {0}. For all /i G R \ {0} we define first 
order and first order symmetric differences by 

\ (Th.x - I) and Sx = S'^ := ^ iSh,x + S-h^x) = ^ (T,,,a - n,-x) , 

respectively, for A G R'^ \ {0} where for all ft. G R we define the shift operator 

Th,\(t>ix) (f){x + hX) 

for functions (p on R*^. We define Sh^ := I and Sq := I. We also adopt the notation 
N = N{-) for a constant N depending only on the parameters given as arguments. 
For basic notions from numerical analysis we refer the reader to |151 120) . 

2. Main results 

To accelerate the rate of convergence of the spatial approximation for a space- 
time finite difference scheme, we will consider the error between the space-time 
approximation and the temporal discretization, the latter of which is a continuous 
problem in the spatial variable. Therefore we begin by considering a discretization 
in time for (11.11) , namely the implicit Euler method. We then replace the differential 
operators with difference operators in this temporal discretization, yielding a space- 
time scheme. We then state our results, the two main results being Theorem 
12.131 which gives sufficient conditions for the existence of the desired asymptotic 
expansion for the error between the space-time approximation and the temporal 
discretization in powers of the spatial mesh size, and Theorem 12.151 which uses 
the expansion to obtain an arbitrarily high order of convergence via Richardson's 
method. 

For each fixed t G (0, 1), we consider 

di 

(2.1) v,ix) = v.^iix) + {LMx) + r ^ + gtii^)) 

p=i 

for a; G 51, X G R'', and i G {l,...,n} with initial condition vq{x) = vq where 
and are second order and first order differential operators given by Ci(j) := 
a^^ {x)DaDii(j) and Mlcj) '■— b°'''{x)Da(t>, for p G {l,...,di}, where the summa- 
tion convention is used with respect to the repeated indices a, /3 G {0,1,..., d}. 
We assume that the given af^ := af^{x) and bf := {bf{xj)'I^Li are real-valued 
and R'^i -valued, respectively, J^i ^ ,B(R'')-measurable functions for w G 57 and 
J G {0, . . . , n} for all a, /? G {0, . . . , d}. The free terms fi := fi{x) and g'^ := .gf (x), 
for p G {1, . . . , di}, are J^i (8) S(i?'^)-measurable functions for every oj G ri, x G R'', 
and i G {0, . . . ,n}. The discretization (j2.ip represents an implicit Euler method 
for approximating the solution to (II. ip in time. Solutions to (j2.1[) with appropri- 
ate initial conditions are understood as sequences of VF2^-valued random variables 
satisfying (|2.1I) in a weak sense in W^j"^- 

As discussed in the introduction, we consider the following degenerate stochastic 
parabolicity condition, necessary for the well-posedness of (|l.ip and hence (j2.ip . 
Note that this is a weaker condition than the strong stochastic parabolicity condition 
which assumes the uniform nondegeneracy of the quadratic form (c/. Assumption 
2.2 in [8 for example). 
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Assumption 2.1. For all oj E fl, i E {1, . . . , n}, x E H'^, and z = (zi, . . . , Zd) E R'^ 

a, 13=1 

that is, the quadratic form is nonnegative definite (positive semi-definite). 

To formulate existence and uniqueness results, as well as estimates, for the solu- 
tion to (|2.1[) we also require some smoothness assumptions on the coefRcients, the 
free terms, and the initial conditions. Let to > be an integer. 

Assumption 2.2. For each u! E ft and i E {0,...,n}, the functions a"''^ and 
the functions a^", af°, and a™ are, respectively, (m + 1) V 2 times and m + 1 
times continuously differentiahle in x for a, f3 E {!,..., d}. For each uj E ^ and 
i E {0, . . . ,n}, the functions bf are to + 2 times continuously differentiahle in x for 
a E {0, . . . ,d}. Further, there exist constants Kj, for j E {0, . . . , m + 2}, such that 

D'af'^ < Kj for j < (m + 1) V 2, 

\D^af\ + 1-0%^"! + \D^a^^\ < Kj for j<m + l, and 

\D^b°'\ + \D^b°\ < Kj for j <m + 2 

for all a, l3 E {1, ■ ■ ■ ,d}. 

For integer / > 0, we define the norm 

n 

:=E^r||0,||f 



1^ < oo. We use the shorthand notation 



and let W2(r) be the space of W^j-valued J^i-measurable processes (f> such that 

1 notation 

p=i 

for functions g = (.g'')pLi- 

Assumption 2.3. The initial condition vo E L^{i},J^o,W™~^'^), the space of J^q- 
measurable W^r^^"^ -valued square integrable functions on fl. The free terms f and 
gP , for p E {1, . . . ,di}, take values in 'W2''~^^{t). Moreover, 

(2.2) JCl := rE ||t;o||',+2 + Ifll+i + MI+i < ^- 

Remark 2.4. For to > d/2, we can find a continuous function of x which is equal to 
vq almost everywhere for almost all € f2, by Sobolev's embedding of C C&, 
the space of bounded continuous functions. Similarly, for each lo E and i E 
{0, ...,«} there exist continuous functions of x which coincide with fi and ijf , for 
p E {!,..., di}, for almost every x E R''. Thus, if Assumption 12.31 holds with 
m > d/2 we assume that vq, fi, and g^ , for p E {I, . . . , di} , are continuous in x for 
alH e {0, . . . , n}. 

For the time scheme (j2.1[) we give the following solvability theorem along with 
an estimate. The proof is provided after some preliminaries are presented in the 
next section. 

Theorem 2.5. If Assumvtions [2A[ \2.S\. and \2.3\ hold, then (j2.1l) admits a unique 
-valued Ti-measurable solution v. Moreover, 

(2.3) Emax||^;,||^<7V/C^ 
holds for a constant N — N{d^ di , m, T, Kq^ . . . , K^n^2)- 
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Now we wish to approximate (j2.ip in space by replacing the differential operators 
with difference operators. Together with (j2.ip we consider, for a finite subset A C 
R"* containing the origin, 

(2.4) vlix) = vt,{x) + (L,Sf (x) + Mx)) r + (Mf^^ti(^) + 5f-i(^)) Cf 

for uj G il, X £ H'^, and i £ {1, . . . ,n} with initial conditions Vq{x) = vq. For each 
? e {0, . . . , n}, the and M^'*''' are given by 

A,/ieA AsAo 

and 

AgA 

for /9 € For all A,^ € A, we assume the given a^^ := a^^{x), := 

and := q^ix) are real-valued and the = are R'^i-valued 

j^i (g) ;8(R'')-measurable functions for every oj £ Q, x £ R'', and i £ {0, . . . , n}. 

In order for to approximate the solution of (|2.ip in space we require the 
following consistency condition, ensuring the difference operators converge to the 
differential operators. 

Assumption 2.6. For all a, /3 £ {1, . . . ,d} and p £ {1, . . . , di}, 

AeAo A,/iGAo 

and 

AeAo /ieAo AeAo ^eAo 

for i £ {0,1,.. .,n}. 

We also place the following additional assumptions on the coefficients of the 
difference operators. 

Assumption 2.7. For all oj £ fl, x £ TV^, and i £ {0, . . . ,n}: 

(i) the functions p'^ > and > for all A £ Aq; 

(ii) for integer c?2 ^ 1 and A G Aq there exist Ti ® B(R.^)-measurable real 
valued functions such that 

(2.5) af'^:=2af^-bf''br = Ea^ar 
for all A, /i G Aq. 

Assumption 2.8. Let I > 1 be an integer. For all lj £ il, i £ {0, . . . ,n}, A £ 
Aq, and k £ {1, . . . ,d2}, the functions and are I + 2 times continuously 
differentiable in x; the functions a^^ are I + 1 times continuously differentiable in 
x; and the functions a^^ , a^*^, a^°, p^, and q^ are I times continuously differentiable 
in X. Further, there exist constants Kj, for j £ {0, ...,/ + 2}, such that 

\D^b^\ + \D^b°\ < Kj for j <l + 2, 

\D^cr^'' I < Kj for j <l + l, and 

jD^af I + iD^a'i^l + \D^a°°\ + \D^p^\ + \D^q^ \ < Kj for j < I 
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for all oj £ fl, X £ R"*, i G {0, . . . , n}, A G Aq, and fee {1, . . . , ^2}. 
Remark 2.9. It is clear that (|2.5p implies that 

XI af^ZAZ^ > 

A,pGAo 

for o; e J7, a; e R'', i G {0, ...,£?}, and za G R for A G Aq. This observation, 
together with Assumption 12.61 implies Assumption 12. II 

Solutions to (12. 4p are understood as sequences of random fields taking values in 
£^(Gh), the space of square summable functions on the grid points Gh, satisfying 
(|2.4p with an £^(G;i)-valued initial condition. The following is a well known result 
which we include for the sake of completeness. Note that by Assumption 12.31 the 
vo, f and gP, for p G {1, . . . , di}, are £^(G;i)-valued processes when restricted to 
the grid G^. 

Theorem 2.10. // Assumptions and \2.8\ hold, then (|2.4p admits a unique 
£'^{Gh) -valued solution for sufficiently small t. 

Proof. The proof of this solvability result relies on the invertibility of (/ — tL'I), 
for each i G {0, . . . , n}, in £'^{Gh) for sufficiently small r. Rewriting the scheme as 
a recursion and using this fact one can construct a unique solution to the scheme 
iteratively. Full details can be found, for example, in [5]. □ 

We observe, however, that (j2.4p is well defined not only at the points of the grid 
but for the whole space. Therefore, we consider (|2.4p on R** and seek solutions 
that are sequences of L^-valued functions. Hence we will use the normal machinery 
from analysis to obtain estimates in appropriate Sobolev spaces for solutions to 
the space-time scheme. Then we will obtain continuous versions of these solutions, 
by Sobolev's embedding, and show that these solutions agree with the "natural" 
solutions at the grid points. 

To aid in achieving this goal one has the following lemma regarding the embed- 
ding W2 C £'^{Gh), the proof of which can be found, for example, in [5]. Recall 
by Sobolev's embedding of W2 into Cb, for I > d/2 there exists a linear operator 
X : W2 Cb such that I(f>{x) = cj){x) for almost every x G R'* and 

sup \I^{x)\<Nmi 
for aU G where N = N{d). 

Lemma 2.11. For all G > "^/^ <ind h G (0, 1), then 

J2 mx)fh''<Nm^ 

xeGh 

for a constant N — N{d). 

With these preliminary considerations in mind, we turn to the main pursuit of 
this paper. To accelerate the rate of convergence of the spatial approximation to 
an arbitrarily high order via Richardson's method we must first prove the existence 
of an asymptotic expansion in powers of the discretization parameter h for the 
error between the space-time approximation and the temporal discretization. Thus 
we prove that for an integer fc > there exists random fields vf^\x), . . . , v^''\x) 
that are independent of h and satisfy certain properties for all i G {0, . . . , n} and 
X G Gh- Namely, that w^"^ is the solution to (|2.ip with initial condition vq and for 
nonzero h, 

(2.6) ^'^^^)=j2Kvl'\x)+Rl''ix) 
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holds almost surely for all x £ Gh and all i e {0, . . . ,n} where v'^ is the solution 
to (|2.4p with initial condition vq and i?'' is an £^(Gft)-valued adapted process such 
that 

(2.7) Emax sup Rj'^{x) ' < N h^'-''+^^ K-l, 

for a constant N independent of h and t. 

We include the following additional assumption on the coefficients of the differ- 
ence operators because at certain points in the proofs to come we will require less 
regularity than is guaranteed by Assumption 12.81 

Assumption 2.12. Let m> be a fixed integer. For A, /i G A, the spatial deriva- 
tives of a^^ and exist up to order (m— 4) VO and for A G Aq the spatial derivatives 
of and exist up to order (m — 2) V and the coefficients together with their 
derivatives are bounded by constants Cm for all u £ fl, x £ R'^, and i G {0, . . . ,n}. 

Theorem 2.13. If Assumvtion [^751 holds with integer I > d/2 and Assumptions 
[M [M m [M [13 andWM hold with 

(2.8) m^m>'ik + A + l 

for integer k > 0, then expansion (|2.6p and estimate (|2.7p hold for h > with a 
constant N ~ N{d, di, c?2, "m, I, T, Kq, . . . , K„i+2, Kq, . . . , Cm, A). //, in addi- 

tion, p'^ = q'^ = for A G Aq, then (|2.6p and (|2.7p hold for all nonzero h. In this 
case, the vanish for odd j < k and, hence, if k is odd, then (j2.8p can be replaced 
with m — m>3k + l + l. 

This theorem follows from the next result, which will also allow us to provide 
higher order estimates for derivatives of the solutions. Taking differences of (|2.6p 
yields 

k 

Sh^xv^x) + J2 ^Sh,xv'i'\x) + 5h^xRI'\x) 



X 



3- 

for any A := (Ai, . . . , Ap) G A^', for integer p > 0, where A" :— {0} and Su 
<5/t,Ai X • • • X 5h,\ ■ Although the estimate for Sh^xRj' '(x) is not obvious, we have 
the following generalization of Theorem 12. 131 

Theorem 2.14. Let the assumptions of Theorem \2.1S\ hold with 

(2.9) m = m>p + 3fc + 4 + / 

for integers I > d/2, p > 0, and fc > with A G A^. Then for h > expansion 
and 

2 



Emax sup 



hold for a constant N — N{p, d, di, d2, m, I, T, Kq, . . . , Km+2, Kq, . . . , Cm, A). 

If, in addition, p'*' = = for A G Aq, then the terms u'^-* vanish for odd j < k 
and, therefore, if k is odd, then (12.91) can be replaced with m > p + 3fc + 1 + / . 

This theorem and Theorem 12 . 131 follow from a more general result that is proven 
in Section H] after some preliminaries are presented in Section [31 Presently we 
formulate our acceleration result, which says the rate of convergence of the spatial 
approximation can be accelerated to an arbitrarily high order by taking suitable 
weighted averages of the approximation at different mesh sizes. 

Fix an integer fc > and let 

k k 

(2.10) f;'^ and i'^ ■.= Y^py-''- 

j=o j=0 
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where solves, with 2^^h in place of h, the space-time scheme (|2.4p with 

initial condition vq. Here (3 is given by (^O: Pi, ■ ■ ■ , Pk) '■= (1,0,..., 0)V~^ where 
is the inverse of the Vandermonde matrix with entries V^^ = 2 '^^^ ' 
forz,j e {!,..., + Similarly, ^ is given by 0oJu---Jk) := (1, 0, . . . , 0)t/-i 
where V^^ s the inverse of the Vandermonde matrix with entries V^^ — 4:^(^^^^'>(i^^) 
for i,j e {1, . . . , fc + 1} where k :— [U . Here [cj denotes the integer part of c. 
Recall that t;'") is the solution to (12. ip with initial condition vq. 

Theorem 2.15. Let the assumptions of Theorem \2.1S\ hold with 
(2.11) m = m>3fc + 4 + ? 

for integers I > d/2 and k > 0. Then 



(2.12) E max sup v1{x)-v'f\x) 



2 



holds for h > with N — N{d, di, d2, m, I, T, Kq, . . . , Km+2, Kq, . . . , Ki^2, Cm-, A). 
If, in addition, = = for A G Aq, then 

(2.13) Emax sup v'^ix) - vf'ix) ' < N\h\^^''+^^ ICf^ 



xeGh 



holds for nonzero h. Moreover, if k is odd, then we only require m = m > 3k + 1 + 1 
in place of (|2.1ip . 

Proof. By Theorem 12. 131 we have the expansion 

^ (0) ^ y Jf^^{2i) + 



2iWj 

i=l 



for each j e {0, 1, . . . , fc} where ''^ := h-^'^+^^R^ Then for r^ Y!1=o 

\3=0 J 3=0 i=l 



,(0) . , 

2i\ ^ 



since X]j=o f^i ~ ^ ^^"^ S^=o /^i^ ^ ^"-"^ each z G {1, 2, . . . , fc} by the definition 
of /3. Now using the bound on W'^ from Theorem 12.131 together with this last 
calculation yields (j2.13p . The result for (|2.12p is obtained in an almost identical 
way and therefore we omit the proof. □ 

Remark 2.16. Note that without the acceleration, that is, when fc = and fc = 1 
in (|2.12l) and (j2.13p . respectively, we have that 

Emax sup Wlix) - Ui(a;)|^ < Nh?K}^^ 



and if = = for A G Aq, then we have 



Emax sup Wl{x) - Ui(a;)|^ < Nh'^K?^ 



in the theorem above. Moreover, these estimates are sharp; see Remark 2.21 in [3] 
on finite difference approximations for deterministic parabolic partial differential 
equations. 
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One can also construct rapidly converging approximations for the derivatives of 
v^^^ by taking suitable weighted averaged of finite differences of w'*. 

Theorem 2.17. Let p > be an integer and let = q'^ = for A G Aq. // the 

assumptions of Theorem \2.13\ hold with 

m^m>p + 3k + 4: + l, 

for integers I > d/2, k > 0, and p > 0, then for A S equation (|2.13|) holds with 
Sh,\v'^ o-nd 6h,\v'^'^^ in place of v'^ and ti^''^ respectively. 

Proof. This assertion follows from Theorem 12.141 in exactly the same way that 
Theorem mis] follows from □ 

We end this section with two examples of ways to choose appropriate a, b, p, q 
and A. 

Example 1. Let A = {eo, ei, . . . , e^} where eo = and is the ith basis vector, that 
is, A is the basis vectors in R'* together with the origin. Then for i G {0, 1, . . . , n}, 
set 

"—a^P and b^"" := 6/ 
for each a, /? e {0, 1, . . . , d} and 

p^ = q^ := 

for each a € {1, . . . Then each spatial derivative Da in (|2.ip is approximated 
by the symmetric difference 5^^ . 

Example 2. Let A again be the basis vectors in R'' together with the origin. For 
J e {0, . . . , n}, set 

a™:=a"" and c^^^' -.^ af 
for each a, /3 G {1, . . . , d} and 

for a S {0, . . . ,d}. We also take Fi ® S(R'')-measurable functions p'^° and cf°' for 
a e {1, . . . , d} such that 

for a e {1, . . . , d}. 

In the next section we make observations that will be used in the proofs of 
Theorems 12.141 and 12.131 which are given in Section S) 

3. Auxiliary Results 

We begin by delivering a proof for Theorem l2.5l For integer to > 0, recall Lemma 
2.1 from [32] taking the parameter p in the Lemma to be p = 2. This Lemma holds 
for all t e [0, T] so in particular we have it for each ir for i E {0, . . . , n}. 

Lemma 3.1. Let^eVF™^^- If Assumvtions \2A\ and WM hold for all multiindices 
7 such that I7I < m, then 

Qlicp) / 2 {D-^cp) D-'C^c^ + V \D-'M'jf dx < N ml 

for a constant N — N{d, di, to, Kq, . . . , Km). 

We use Lemma [01 to obtain estimate (|2.3p . The existence of a solution to (|2.ip 
will follow from the vanishing viscosity method. 
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Proof of Theorem \2.5\ We first assume that a sufficiently smooth, solution to (|2.ip 
exists and obtain estimate (|2.3p for a constant N independent of r. We begin by 
obtaining an expression for the square of the norm for the solution to the time 
scheme. Then we estimate the supremum of the expectation of the square of the 
norm and in particular we show that this quantity is finite. With these observations 
in place we are then able to estimate the expectation of the supremum of the square 
of the norm. 

For a multiindex I7I < m, considering the equality o? = 2a{a — b) — \a — 5p 
we note that (j2.4p implies 



\D-'v.,\\l ~ \\D^v,^i\\l = 2 iD^v,,D^ {v, - v.,^i)) - \\D^ {v, - 



2 [p-^v^.D-^ {C,v, + /,)) T - [v, - 

di 



p=i 



2 {D-^v^.D-^ {Uv, + /,)) r - (At-, + /0||o 



2 Ji 



d ^ 

Yd-^ (A^r-i-.-i+sf-oef 
p=i 

Summing up over i from 1 to j and over I7I < to, we have 
(3.1) 
where 



\\v,\\i<\\vo\\l + H,+I, + J,, 



i=i 



j di 



and 



i=i p=i 



Js = E 



Yd-- (A^ri«.-i + 5f_i)ef 
p=l 

By an application of Ito's formula, for each tt, p e {1, . . . , di} one has that for 
all i e {0, . . . , n - 1} 



ciViC^Vi = «i - - -YZ- Yr + Tx.p 



for 



where k(s) is the piecewise defined function taking value k(s) = i for s G [ir, (i+l)r) 
and where x-np = 1 when tt = p and otherwise. Thus we can write Jj — jj^' + jj^' 
where 



7(1) 



E 



E^™(>'f-i--i+5f-i] 

p=l 
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and 

jm r j2 (i?"(A^:j(.)«.(.)+5:(.)),i?'"(K(3)X'«(s)+5:(.)))dy-''(.). 

Now observe that, for each i G {1, . . . ,n}, by Lemma [3. II we have 
H,+//^<NrJ2\\DaVo\\l + Nrj2 E (Q7(«.) + PVJo + ll^"3-illo) 

a=0 i=l |7|<m 

(3.2) < Nr WvoWl^, + Nrj2{\Ml + WMI + h^-iWl) , 

1=1 

where N — N{d, di,m, Kq, . . . , Km+i)- Note that we only require 6" to have 
bounded derivatives of at most order m + 1; the other coefficients only need to 
have bounded derivatives of at most order m at this stage. Here the initial con- 
dition vo enters, estimated in the W™'''^-norm, due to the displacement caused by 
the discretization in time when we consider the quadratic form from Lemma 
13.11 Thus inequality p.ip becomes 

j 

(3.3) \\v,\\l < NT \\vo\\l+, + A^r £ + + \\g,.^\\l) + I, + jf . 

1=1 

Since Vi, M'^Vi, and are all J-^-measurable and ^^^-^ is independent of Ti for 
I € {0, . . . , n}, we have that 

= 2^ {(^'^^-1. {Ml,u,-i + g")) E (^f I J-,_i)} = 0. 

i=i p=i 

(2) 

Similarly, we see that EJj ' = since the expectation of the stochastic integral 
is zero. Therefore, taking the expectation of p.3p and the sum of / and g over 
i € {0, . . . , n}, we have that 

j 

(3.4) E \\v,\\l < N (rE \\vo\t^, + ml + Wl) + NtEJ^MI 

1=1 

for each j G {1, . . . ,n}. Applying a discrete Gronwall lemma to (I3.4p we have 

E \\v,\\l < N (rE WvoWl^, + ml + Ml) (1 - Nry^ 

and, since 

(l-«.r'^(l-A'^)"<(l-A'|)""<Ce-, 

we have the following estimate for the supremum of the expectation of the square 
of the norm 

(3.5) maxE \\v£^ < N (rE \\vo\\l+, + ml + Ml) 

for a constant N = N{d, di, m, T, Kq, . . . , Km+i)- In particular, we can use p.Sp to 
eliminate the last term on the right-hand side of p.4p by bounding it with terms 
already appearing on the right-hand side p.4p . 
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Next we approach the estimate for the expectation of the supremum by first ob- 
serving how to bound the I and J^^^ terms appearing in p.3p using the Burkholder- 
Davis-Gundy inequahty. For J*-^' we have 



(2) 



E max I J, 



<cy: e. 

7T,p—l 



1/2 



y / y ■> rn m 

. 1/2 



TT^p—l 



< 



C y EmaxV7|l7Wfwi 

^ — ^ i<n 



-5flL 



7T,p—l 
T 



1/2 



ds 



where C is a constant independent of the parameters and functions under consid- 
eration and is allowed to change from one instance to the next. Therefore, 



Emax 



J 



(2) 



<rfiC VrEmaxllTWfwi 



(3.6) 



p=i 



E Erii-M:(.)««(.)+<(.£i^-(.)-<(.) 



ds 



by Young's inequality. We observe that the second term on the right-hand side of 
p.6p can be estimated by 



ds 



< 







2 










dsj 









<iVrE5;]||x;,||^^, + 7V[ 



Im+l 



using the tower property for conditional expectations. Further, the first term on the 
right-hand side of (I3.6p is bounded from above by the sum over all i e {1, . . . , n} 
and can be estimated by the same quantity. Combining these estimates and using 
3)) with m + 1 in place of to we see that E max | J*^^) | is estimated by 



(3.7) 



E max 



(2) 



< 



^(^E||^;o|lL+2 + I/t+l + I.C+l 



for a constant N — N{d, di, m, T, Kq, . . . , K„i+2)- Moving on to /, we note that 



i=i p=i 

= 2^ r du;''(s). 
p=i "^0 
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Applying the Burkholder-Davis-Gundy inequality once again, we obtain 




and then using Young's inequality followed by p. 51) with to + 1 in place of to we see 
that Eniax|/| is estimated by the same quantity appearing on the right and side 
of (1^ . 

Returning to p.3p . taking the maximum followed by the expectation, and using 
the estimates for the expectation of the suprenium of jJ^^-*! and |/|, we see that 

(3.8) Emax < N (rE \\vo\\l+, + ml+, + MI+,) = NK^, 

holds with a constant N = N{d^ c?i, m, T, i^o, . . . , Km+2), thus establishing (|2.3p . 
Next, we use the vanishing viscosity method to show that (|2.ip admits a solution. 

For e > 0, we let £f :— Ci4> + eA0 where A :— X)q=i DaDa is the Laplacian. 
Notice the leading coefficient of the operator £f is given by a"^ :— a^^ + exaji, 
where Xap = 1 for a = /3 and zero otherwise. We then consider the equation 

di 

(3.9) vt = vU + (>cf < + f.)T + Y, {MUvt-i + gU) Cf 

p=l 

for each i S {1, . . . ,n} with initial condition Wq — vq. Proving the solvability of 
p.9p reduces to solving, for each a; S fi, the elliptic problem 

(3.10) Avt = 
for each i € {1, . . . , n} with free term 

di 

F.:=vU+rf, + J2^' i^li^U + all) 
p=i 

and operator 

A := {I-rCD 
where / is the identity. That is, we claim that Ai is 

(i) bounded, i.e. \\Ai(t)\\^ < J'^||0||m+2 for a constant K, 

(ii) and coercive for sufficiently smafi t, i.e. {Ai<f),(f)) > ||!<?!'|!m+2, 

for all (j) e W^+'^ for every i G {1, . . . , n}, where (•, •) denotes the duality pairing 
between W™'^'^ and based on the inner product in W^~^^ . We will obtain the 
existence of a solution vf to (|3.10p for each i G {1, . . . ,n} via Galerkin approxima- 
tions (of course, in this instance, one could also use the Lax~Milgram Theorem). 

For integer p > 0, let Ep be the p-dimensional subspace of W™^^ spanned by 
the first p elements of {ej;j G N}, a collection of vectors from W™^^ forming an 
orthonormal basis for W™^^- We seek an approximate solution (p^ e Ep to 

for each k G {1, . . . ,p}. Rewriting (p^ — cj,ej for coefficients c^, where the summa- 
tion convention is used with respect to the repeated index j G {1, . . . ,p}, we see 
that is an approximate solution if and only if is W^'^^ -valued and satisfies 
the system of ordinary differential equations 

c^p {Aiej,ek) = {F.i,ek) 
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for each p. We derive the estimate 

From this we see that {AiCj, Ck) is invertible and thus a solution c^, and hence an 
approximate solution 0^, exists and moreover we have the estimate 

E||(/)f ||,„+i < 2e-iE||^,||,„ 

uniformly in p. Thus, there exists a e W^^"^ and a subsequence pk such that 
1^^'° — weakly in W™"*"^. Therefore, for each i £ {1, . . . ,n} there exists a vf 
satisfying (I3.10p and, moreover, this solution is easily seen to be unique. Hence, we 
construct a unique solution to (I3.9P iteratively. 

Using the existence and uniqueness to the elliptic problem in each interval, we 
note that 

di 

vo + T/i + ^ {M'o^o + g'o) e 
p=l 

by the assumptions on vq, /, and g^, for pG {l,...,c?i}, and therefore there exists 
Avf e W^^'^ satisfying 

A^v\ = vo + r/i + ^ (Mgt'o + g'o) ei- 

p=0 

Further, assuming that there exists a e W^2"^^ satisfying (EH), we have that 

dl 

by the induction hypothesis and Assumption [531 and therefore there exists a wf_|_]^ S 
^2™"^^ satisfying ([221). Hence, we obtain = (vf)"^^ such that each vf G 1^2"""^^ 
satisfying (j3.9p . 

Finally, we observe that the estimate p.Sp can be obtained for the solution to 
p.9p in a similar manner. In particular, this gives a uniform estimate in e for the 
solution to p.9p . Therefore, there exists a subsequence — > and a WJ"-valued 
J^i-measurable Vi such that wj^*" converges weakly to as fc — > 00 in W™ for each 
i G {1, ■ . ■ This V — (t^i)"^! is the solution to (|2.ip and is easily seen to be 
unique. □ 

Remark 3.2. We consider an implicit scheme where the operators £ and take 
values at the points of the time grid. The displacement observed in p.2p caused by 
the discretization in time can be avoided if we consider a modified implicit scheme. 
Namely, if we consider operators defined to be the average over the intervals defined 
by consecutive points of the time grid, as in [7], we could then take M^vq ■= 0. 
However, we believe such a scheme would be less practical from a computational 
standpoint. 

The lemma below is given in [1] for all t G [0, T] so, in particular, we have the 
following for each ir for i G {0, . . . ,n}. This lemma plays the role of Lemma [3.11 
for obtaining the estimate for the space-time scheme. 

Lemma 3.3. Let G Wj^^. // Assumptions lS77\ and [KB\ hold, then for all multi- 
indices 7 such that I7I < I, then 

d 



Q7(</.) = / 2 (Z?^</>(x)) D'<L':<t>{x) + li^^Mf 



\x<N\\cj>\\l 



for all i G {1, . . . , n} for a constant N = N{d, di, d2, 1, Kq, . . . , Ki^i, A). 
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We are now able to give an estimate for solutions to the space-time scheme that 
is independent of h. Recall that for integer Z > 0, we define the norm 



and let W2(t) be the space of W^j-valued J^i-measurable processes (f> such that 
Ufi < oo. 

Theorem 3.4. // Assumptions \2.3l \2. 7[ and \2.8\ hold, then for each nonzero h 
equation (|2.4p admits a unique W^-valued J-i-measurahle solution . Moreover, 
satisfies 

,2 



(3.11) 



Emaxllvfll" < NlCf 



for a constant N — N{d, di, d2, 1, T, Kq, . . . , Ki^2, A). //, in addition, p 







for A G Aq, then (j3.11l) holds for all nonzero h. 

Proof. That (j2.4p admits a unique L^-valued solution follows immediately from the 
considerations in the proof of Theorem 12.101 In particular p.lip can be achieved 
easily for a constant N depending on /i, so we see that the solution v'^ is "Valued 
and J^i-measurable. To achieve p. lip for a constant N independent of h (and r) 
follows almost immediately from the derivation of the estimate (|2.3p in the proof 
of Theorem 12. 5[ using Lemma 13.31 in place of Lemma 13.11 
We obtain 

(3.12) Wvff^ < \\vof, +Hj+Ij+Jj, 

in the same manner as p.ip . with 



i=l 



j d 

El 

i=i p=i 



9li 



and 



i=l 



Then by an application of Ito's formula, we rewrite Jj — J-^' + J)^' using 



7(2) 



2 



p=i 



and 



J. 



(2) 



TT.p—l 



where yP{t), for t e [Q,T], is defined in the proof of Theorem (^3]) . 
Now observe that, by Lemma [3?3l for each i G {0, . . . , n}, we have 

H, + Jf^ < r J2 \\Sk,xVo\\' + Nrj2 (QZK") + {D'' , D'' f.) + \\D' 

AGA i=l 

< NT ||«o||f+i + iVr^ (}\vlt + 11/^11? + 



9i-i\ 



i=l 
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for /i > where N = N{1, d, d2, Kq, . . . , Ki^i,A). Again, here the initial condition 
?;o enters, estimated in the H/^^^-norm, due to the displacement caused by the 
discretization in time when we consider the quadratic form Q'^ from Lemma l3.3l If, 
in addition, p"** = q"** = for A G Aq, then this last calculation holds for all nonzero 
h. Thus, inequality p.l2p becomes 



(3.13) \\vf\\' < NT ll^ollti +NtY, (Ih'll? + \\h\h + ll.9'^-ill?) + + J: 



(2) 
' 



Since E/j = and ^jf'^ = 0, taking the expectation of (|3.13l) and taking the sum 
of / and g over i e {0, . . . , n}, we have that 



(3.14) E\\v^l <N[rE\\vo\\l, + ml + Ml)+NTEj2\Hl 

i=i 

for each j e {!,...,«}. Applying a discrete Gronwall lemma to (I3.14p we have 
E llz;^"!!^ < N (rE \\vo\\l, + [/if + M?) (1 - Nt^ 



and thus 



(3.15) maxE II' < TV (rE \\vo\\^^, + [/I 

i<n ^ \ 



I + Ml 



for a constant N — N{d, di,d2,l,T, Kq, . . . , Ki^i, A). In particular, we can use 
p.lSp to eliminate the last term on the right-hand side of p.l4p by bounding it 
with terms already appearing on the right-hand side p.l4p . 

The terms / and J*-^-* are estimated as in the proof of Theorem 12.51 with M'^'P 
in place of Ai^, using the Burkholder-Davis-Gundy inequality. In particular, we 
have 



E max \ Ii\ -t- E max 

i<n i<n 



(^Ellvoll 



<iV(rE||z;nl|;+2 + I/lf+i + M'+i 



for a constant N = N{d, di, ^2, ^ T, Kq, . • . , A). Thus, returning to (|3.13p and 
taking the maximum followed by the expectation we have that 

Emax \\v1\\' < N (rE ||«o||?+2 + I/l'+i + Ifflf+i) = NJCf < 00, 

for a constant N ~ N{d, di,d2,l, T, Kq, . . . , K1+2, A). □ 

For the convenience of the reader we record the following lemma, found in [3]. 

Lemma 3.5. Let (f> € M/^f^^ and xjj e W^2^^ for an integer p > and let X, ^ G Aq. 
Set 

dx(j) := XWj4> and := dxd^. 

Then we have 

(3.16) ——Sh,xH^) = e^di+'^ix + hex) ae, 



(3-17) 77n:T;7'5A</>(a;) I epdl+^(p{x + hex) de 



QP 1 /•! 

and 
(3.18) 

5x5^4,{x) f {9,dx - e2d^)Pdx^ij{x + h{eA - 02^l)) de^ de2 



{dh)p '^-^ ' 4 
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for almost all x G R"^ for each h E R. Thus 

QP 



(3.19) 
and 

(3.20) 

where 
(3.21) Bp 



{dh)P 



QP 



B. 



.0 



5P 

{dh)P 



h—0 r—0 



if p is odd 

1 if p is even 



A — 



if p or r is odd 

^ (r+iy.[p-r+iy. '^fP S^S" 



Furthermore, for integer I >0, (f) (£ Wl"^^^', and € one has 



(3.22) 



J=0 



+ 



< 



p+1 



(p + 2)! 



(3.23) 

and 
(3.24) 



+1, 



< 



\h\P+^ 
(p + 2)! 



j=0 r=0 



l+p+3 ' 



w/iere TV = iV(|A|, \p.\,d,p). 



For integers ^ > and r > 1, denote by the Hilbert space of functions (p on 
R"* such that 



(3.25) 



and set W, 



r,h — E X • • • X 5h^xA\?i < ^ 

Xl,...,\r-GA 



h,2 



Remark 3.6. Formula p.l6p with p = and Minkowski's integral inequality imply 
that 

\M\\o<\\9x<t>\\o- 

By applying this inequality to finite differences of (f> and using induction we can 
conclude that Wj^"" C VF^'j- Further, for any e Wj^*" we have 



where N depends only on |Aop = J2 



AgAq 



P and r. 



We now use the preceding observations to obtain estimates in appropriate Sobolev 
spaces for a system of time discretized equations. Here we use the summation con- 
vention with respect to the repeated indices A, /i € Aq. For i e {0, . . . , n}, let 

4°' ■.= 4''d,d,+ {p^^q^)d,, 
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for each {1, . . . , di}, and for an integer p > 1 let 
p 



-j+i 



Br, 



AO 



,0A\ 



M 



(p)p 



B 



O. 



Hp) ._ 



P U^PfiP+'^ 

h3 



and 



hi 



U)p 



J'- 



where Ap,j and Bp are defined by p.2ip . By AssumDtfon l2.61 we have that 

f^P = TWf for all i G {0, For p > 1, the values of c'fU and 

are obtained by formally taking the pth derivatives in h of L^cj) and Ml^'^cj), 



and M 



M 

respectively, at /i = 0. 

Remark 3.7. Let Z and p be nonnegative integers. For G W^''''^'''' and -0 G W^'''''^'''' , 
under Assumptions 12.21 and l2.12] with m = m we have that for I < to, 



7^: 



h{p)p. 



(3.26) 
and 
(3.27) 

for each p e {1, . . . , by for a constant TV = iV(p, d, to, iiT/, C„, A) 

For integer fc > 1, the sequences of random fields v^^\ . . . , w*^*^^ needed in 



will be the embedding of random variables taking values in certain Sobolev spaces 
obtained as solutions to a system of time discretized SPDEs. Namely, as the solution 
to 



,(p) 



(p) 

'4-1 



0"),,(p-j) 



(3.28) 



for p G {1, . . . , fc} where ~ p{p — 1) . . . {p — j + l)/j\ is the binomial coefficient 
and v^-^^ is the solution to (|2.ip with initial condition vq. 

Theorem 3.8. // Assumptions \EJl [13 and WI^ hold with m = to > 3fc, 
then the system (j3.28p admits a unique solution i/^^-* , . . . , v^^'^ for initial condition 
— such that is ^P -valued Ti-measurahle. Moreover, for 



,(1) 



each p £ {1, . . . , fc} if v^^^ is a solution, then 

2 



(3.29) 



E max 



,(P) 



< NKi 



m—3p 



holds for h > with a constant N — N{d, di,m, k,T, Kq, . . . , Km+2,Cm)- If in 
addition, p'^ q-^ = for A £ Aq, then p.29p holds for all nonzero h and — 
for odd p < k. 



20 



E.J. HALL 



Proof. For convenience let 



and 



where we will write G^^) = Y^tU G^^^''- Observe that for each p £ {l,...,k} the 
equation for i^(P) in p.28p depends only on v*^^' for j < p and does not involve 
any of the unknown processes v^^'^ with indices j > p. Therefore, we shall prove 
the solvability of the system and the desired properties on i^^^-* recursively using 
Theorem O 

By Theorem 12. 5[ is W™-valued, J^^-measurable, and satisfies p.29p with 
p = 0. Observe that 



< TV 



,(0) 



for a constant N = N{m, Cm), recalling the constant Cm from Assumption 12 . 1 2 1 in 
this instance since we only require that (m — 2) V derivatives of the coefficients 
exist and are bounded, and further that M^^vf^ = 0. Therefore, by Theorem l2.51 
there exists a unique M^™~^-valued J^^-measurable v^^^ satisfying (j3.28p with zero 
initial condition. Moreover, the estimate p.29p is clearly satisfied in for p = 1. 

Now we induct on p, assuming that for to > 3fc > 2 and p G {2, . . . , fc} we have 
unique solutions v^^^ . . . , ^/^^'^^^ satisfying the desired properties. Observe that 





< N 


(p-i) 






m — 3p+l 







< N 



m — 3(p— 1) 



m-2(p-l) 
(1)P,,(P-1) 



for a constant N = N{m,Cm), by Assumption 12.121 and that " v, 

iv-j) 



= 0. 



Further for j > 2, if j is even, then 





< N 






2 Z 


m — 3p+l 




m 



< N 



in-3p+2(j"-l) 



m-3{p-j) 



and 



2 



p=i 



m-3p+(2j-l) 



< N\\u 



iP~3) 



m-3{p-j) 



or if j is odd, then 



and M 





< N 




< N 




m — 3p+l 




m-3p+2(j"-l) 



m-3(p-j) 



U)p,,(p-j) 



= 0, for constants N — N{m,Cm), by Assumption 12.12] There- 
fore, by the induction hypothesis, F^p^ is -valued and J^^-measurable, G*^^^ 
is Wj7^~^''^^-valued and J^^-measurable, and 

(3.30) [^(P^r + iG^P^r <NICl 



m — 3p+l 



/■™-3P+V^A c^A nip) d W'"-3P+1( 

there exists a W™~'^^-valued J^i-measurable i/^p-' satisfying (|3.28p with zero initial 



That is, e W""^p+^(t) and G^p) G W™^''^^' (r). Thus by Theorem 



condition. Moreover, Theorem 1 2 . 51 vields the estimate 



E max 


(P) 


^ < iV ( 


pip) 


2 

+ 


G-Cp) 


' ) 






m — 3p \ 




m — 3p+l 




m — 3p+l y 
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for a constant N = N{d, di,m,k,T, Kq, . . . , Km+i,Cm)- Combining this with 
p.30|) yields (I3.29|) . We remark further that the uniqueness of each z/^^^ follows 
from the uniqueness of the solutions obtained from Theorem 12.51 

Note that M^p^p = for odd p < fc by (1X1^]) . Assume, in addition, that 
= = for A e Aq. Then also C^p'^ = for odd p < fc by ((3l9)) and dSSO]). 
Therefore, F'^' = and G'-^-' = 0, which implies v^^^^ = 0. Assume that k > 2 and 
that for an odd p < k we have i^^^^ = for all odd j < p. Then C^P'i'liyW) — and 
;V^(p-i)Pj,(j) ^ for all j £{!,..., p} since either j or p-j is odd. Thus F'^p'' = 
and G(P) = 0. Hence = for aU odd p < k. □ 

For integers A: > let i>^^\ . . . , i/^'^^ be the solutions to the system (|3.28p with 
zero initial condition coming from Theorem 13.81 Let 

k 



(3.31) 



,(0) 



0) 



for i g {1, . . . , n}, where v'^ and i^^^-* are the unique VF™-valued solutions to (j2.4p 
and (|2.ip . respectively, with initial condition vq. 

Lemma 3.9. Let v'^''^ be defined as in (IX?T|) . // Assumvtions [KT[{EM fOi [Ol 
and\KT^ hold with m = m>3fc + 4 + ? /or integers k > and I > 0, f/ien 



g — and x^'^ is W2-valued J-i- measurable such that 



Emax 



T,h 



< oo 



T.h 



and r"^'^ satisfies 
(3.32) vl'' - 

for i £ {1, . . . ,n}, where 

and 



j=0 



di 

E 

p=i 



M. 



h,p tJi 



(3.33) 



h.p 

0.-1 



Moreover f G W+^{t) and q'^'P G W'+i(r). 

Proof. First recall that by Theorem 13.41 the solution to the space-time scheme v'^ 
is W^-valued J^i-measurable and satisfies estimate p. Ill) owing to Assumptions 12.31 
12.71 and 12.81 with m = I. By Theorem 12.51 z^'"-* is VK2"'^8'lued J^^-measurable and 
satisfies estimate p.3|) by Assumptions 12. ll 12.21 and 12.31 with m = I. By Theorem 
13.81 the ly'^^^ are W^j'^^^lued J"i-measurable processes satisfying estimate p.29p for 
all j e {1, . . . , fc} owing to Assumptions [2Tl [2?2l [2?3l and [2?T2] with m = m = I. 

A- 

2 



Thus r^''' is M^i-valued J", -measurable and satisfies 



E max 

Kn 



T.h 



< OO 



for m = m = / > 3A: 

One can easily show that r^''* satisfies ((53^ using (gj]), (113), and by 
noting that we can rewrite f'' and g''''' as 



rh ^ j^h^m 
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and 



for 



k ■ k / 



i=0 j=0 



3'- 



1=0 i=l ■' 

EEf^^^^'-^^"^ 

i=i j=o 



k 

EE 
EE 



and 



k—i 



hi 



k i 



E^E^-^^^-^^-'^^^EE- 



i! — ' 7 



._A/f(j)Pj^(«-i) =. J 



where summations over empty sets are taken to be zero. 

To prove the last assertion, we note that, by assumption, m — 3j > I + k — j + 4 
for j G {0, 1, ... , k}. Thus by Lemma 1531 for j g {0, 1, . . . , fc}, i e {0, . . . , n}, and 



and 



^r>h{k-j) (j) 

i i 


< c 










l+k- 




< c 








l+l 




l+k 



< c 



< c 



,0) 



m-3j 



m-3j 



for a constant C independent of the functions and parameters presently under 
consideration that changes from one instance to the next. Therefore f'^ G W2"^^(r) 
and fl'*''' G W^+\r). □ 



In the next section we prove the expansion results. 



4. Proof of Main Results 



We prove a slightly more general result which implies Theorem 12.141 and hence 
Theorem 12.131 Here we suppose that m = m everywhere. 

Theorem 4.1. Letv'^'^ be defined as in (I3.3ip . If Assumvtion [E7S\ holds with integer 
l>0 and Assumptions\K^\EM KM \MM \^ and\EM hold with 

(4.1) m = 3fc + 4 + ; 

for integer k > 0, then, for h > 0, 



(4.2) 



Emax 

i<n 



T,h 



holds with a constant N = N{d, di, d2, m, I, T, Kq, . . . , Km+2, Koi ■ ■ ■ ^ Ki^2, Cm, A). 
//, in addition, p"^ = = for A G Aq, then (|4.2D holds for nonzero h and it 
sujfices to assume m > 3k + 1 + I in place of (14.11) . 



HIGHER ORDER APPROXIMATIONS FOR DEGENERATE EQUATIONS 23 



Proof. By Lemma [231 we have that f e W^+^r) and q'^ p e W?,+^(t). Thus, by 
Lemma 15^ and Theorem 13.41 



E max 


T.h 


;<iv( 


i<n 







(4.3) 



for a constant N — N{d, di,d2,l, T, Kq, 
and by Remark 13. 261 we have 

<N\h\''-'+^ 



i i+i 



, krn+2)- Then by gH) for j£{0,...,k} 



(4.4) 

and 
(4.5) 



< iV|/i| 



k-j+i 



l+k-j+4 



< N\h 



k-j+i 



l+k-]+3 



<N\h\ 



k-j+1 



,0') 



,(j) 



m — 3j 



in — 3 J — 1 



Now using Theorem 13.81 we see that 



which, when taken together with (|4.3p . imphes (j4.2p . 

If, in addition, p'*' = = for A G Aq, then as in Theorem 13.81 it fohows that 
for ah odd j < k. If fc is odd, then clearly i^^*^) = and (|4.4p holds for 
j = k and also for j > k — 1 and we need only m = 3fc + 1 + L We mention further 
that v^^^ — for all odd j < fc in the case = q'^ = for A G Aq also follows from 
(j4.2p . now valid for all nonzero h, since and i^"'' are the L^-valued solutions to 

due to the uniqueness 



(j2.4p with initial condition vq and we must have 
of solutions. 



□ 



We have the following corollary to Theorem 14.11 which implies Theorem 12 . 141 and 
hence Theorem 12.131 Let K^'^ := Ir^''*, where I is the embedding operator from 



[^m and notice that ' £ f{Gh) for aU i e {0, . . . , n}. 

Corollary 4.2. // in Theorem \4. 1\ we have I > p + d/2 for integer p > 0, then, for 
A e AP and h> 0, 

2 



Emax sup 



holds with a constant N — N{d, di, ^2, m, I, T, Kq, . . . , K„i+2, Kq, 



Proof. By Sobolev's embedding of W2 ^ into Cb and Remark 



E max sup 



,Rj'^{x) <7VEmax v^'^ 

i<n 

which, together with Theorem 14. 11 yields the desired result 



■ • , Ki+2,Cm., A). 

we have that 
2 



< A^Emax 

l—p,p,h i<n 



□ 



It only remains to explain how the corollary implies Theorem 12.141 and hence 
Theorem 12.131 Using Sobolev's embedding, we find continuous versions of the L^- 
valued solutions to the space-time scheme, the time scheme, and the system of time 
discretized equations. Then we notice that the restriction of the L^-valued solution 
to (|2.4p onto Gh agrees with the unique P{Gh) valued solution. This argument 
can be found in [5], nevertheless we reproduce it here for the convenience of the 
reader. 

For I : W2 — > Cb from Lemma 12.111 Theorem 12.141 follows by considering the 
embeddings :— Iv^ , where is the unique L^- valued solution to (j2.4p with 
initial condition vq, and w'-'-' :— for j € {0, . . . , fc}, where v'^^^ is the unique 

L^-valued solution to (|2.ip with initial condition vq and the processes z/^^^ , . . . , v^^'' 
are the solutions to the system of time discretized SPDEs p.28p as given in Theorem 
By Theorem l3.41 is W^2"'^9'lued and J^^-measurable for alH S {1, . . . , n}. For 
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each j e {1, . . . , k}, the v^^^ are M^'"~'^''-valued processes by Theorem 12.51 Since 
I > d/2 and m — 3k > d/2, the processes u'* and v^-^^ are well defined and clearly 
p.3ip implies (I2.6P with in place of v'^. That is, we have the expansion for a 
continuous version of the L^-valued solution. 

Next we show that the restriction of the L^-valued solution to the grid Gh, a set 
of Lebesgue measure zero, is indeed equal almost surely to the unique i!'^(G;i)-valued 
solution that one would naturally obtain from (12.41) . That is, we show that 

(4.6) 

almost surely for all i £ {l,...,n} and for all x e where v'^ is the unique 
J^i-measurable €^(G;i)-valued solution of (12. 4p from Theorem 13.41 Therefore, for a 
compactly supported nonnegative smooth function cf) on R'^ with unit integral and 
for a fixed x G Gh we define 

My) ^=<^(^) 

for y G R'' and e > 0. Recall, by Remark 12.41 that we can obtain versions of vq, 
/, and g'', for p e {1, . . . , di}, that are continuous in x. Since v'^ is a L^-valued 
solution of (j2.4p for each e, almost surely 

/ v':iy)My) dy= [ vtAy)My) dy + r [ {iJ^v'i + h){y)Uy) dy 

J'R.'^ JR'' JR<i 

+ Eef/ (Mf_:^i*ti+gri)(y)0e(y)dy 
p=i -^R" 

for each i € {1, . . . ,n}. Letting e — 0, we see that both sides converge for all 
i G {1, . . . , n} and w G fi. Therefore, almost surely 

di 

vlix) = vt,{x) + (L'lv'^ix) + Mx)) r + {Mt^vt,{x) + gl.ix)) 

for all i G {l,...,n}. Moreover, by Lemma [2.111 the restriction of w'', the con- 
tinuous version of i/'*, onto G^ is an ^^(G;i)-valued process. Hence, (j4.6p holds, 
due to the uniqueness of the J^i-measurable ^^(G/i)-valued solution of (|2.4p for any 
Jx)-measurable ^^(G/i)-valued initial data. This finishes the proof of Theorems l2.14l 
and [2131 
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